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^ ■ ABSTRACT 

Context. The observed dynamical mass-to-light [M / L) ratios of globular clusters (GCs) are systematically lower than 
the value expected from 'canonical' simple stellar population models, which do not account for dynamical effects such 
as the preferential loss of low-mass stars due to energy equipartition. It was recently shown that low-mass star depletion 
can qualitatively explain this discrepancy for globular clusters in several galaxies. 
. Aims. To verify whether low-mass star depletion is indeed the driving mechanism behind the M/ L decrease, we aim to 

' predict the M/Ly ratios of individual GCs for which orbital parameters and dynamical V-band mass-to-light ratios 

VpH' M/Lv are known. There is a sample of 24 Galactic GCs for which this is possible. 

Q , Methods. We use the SPACE cluster models, which include dynamical dissolution, low-mass star depletion, stellar evo- 

5_( . lution, stellar remnants and various metallicities. We derive the dissolution timescales due to two-body relaxation and 

disc shocking from the orbital parameters of our GC sample and use these to predict the M / Ly ratios of the individual 
' GCs. To verify our findings, we also predict the slopes of their low-mass stellar mass functions. 

Results. The computed dissolution timescales are in good agreement with earlier empirical studies. The predicted M/Lv 
are in Icr agreement with the observations for 12 out of 24 GCs. The discrepancy for the other GCs probably arises 
because our predictions give global M/ L ratios, while the observations represent extrapolated central values that are 
^ ' different from global ones in case of mass segregation and a long dissolution timescale. GCs in our sample which likely 

CO , have dissimilar global and central M/L ratios can be excluded by imposing limits on the dissolution timescale and 

OO ■ King parameter. For the remaining GCs, the observed and predicted average M/Lv are 78ln% and 78 ± 2% of the 

I canonically expected values, while for the entire sample the values are 741^% and 85 ± 1%. The predicted correlation 

between the slope of the low-mass stellar mass function and M/ Ly drop is found to be qualitatively consistent with 
. observed mass function slopes. 

' Conclusions. The dissolution timescales of Galactic GCs are such that the ~ 20% gap between canonically expected 

' and observed M/Lv ratios is bridged by accounting for the preferential loss of low-mass stars, also when considering 

, individual clusters. It is concluded that the variation of M/L ratio due to dissolution and low-mass star depletion is a 

plausible explanation for the discrepancy between the observed and canonically expected M/L ratios of GCs. 
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1. Introduction stellar populati on (SSP) models that assume a 

^, . . , . , T 1 ..x/.N . r canonica l IMF dR eikuba et al.' '2007*; 'Kruijssej [20081; 

The topic of dynamical mass-to-light (M/L) ratios of KruiissenJ^Lamers 2008; Mieske et al. 200i). 

old compact stellar systems has attracted increasing at- _ a consequence, the dynamical M/L ratios of UCDs 

tention during recent years_tMcLaughhn & van der Mard ^„ ^^.^.^^.^ , -^^^ 1^,^^ ncs. at 



— — ir— ° nu ^'Vs^ -^-i , , TT-i, are on average about twice as large as those of GCs, at 

2005':'Hasegan et al.112005: Rcikuba et al. 2007: Hilker et aP pnrr,n«rpblp ^pf«lhrltlp« 



't^^''': '''' ivci..LLu^ c <.i ^uu , , comparable metalhcities 

2007; E vstigneeva et al. 2007; Da brmghausen et al.l 120081; ^ 
Mieske et al.l l2008nK ruiisscn 20 08|; IBaumgardt fc Mieskel 



2008 HChilingarian et al.ii2008< : iForbes et al.ll2008[) . The out- Regarding GCs, a viable solution to obtain lower M/L 

come of these studies can be summarised as follows: '"^tios^ is a deficit m low-mass stars with respect to a 

canonical initial mass function (IMF) (IKrouDall2001h . This 

— For the mass regime of ultra-compact dwarf galaxies is known to arise naturally from two-body relaxation in 
(UCDs, M > 2 X 10^ M0), dynamical M/L ratios tend star clusters, which causes a depletion of low-mass s tars 
to be some 5 0% above predictions from stellar popu- (IVesperin i fc He ggielllQQTt iB aumg ardt fc MakiMll2003( ). In 
lation models dPabringhausen et aLll2008l : [Mieske et al] iKruiissen (2008) it has been studied how the preferen- 
I2OO8I : iForbes et al.ll2008f ). tial loss of low-mass stars due to dynamical evolution of 

— For the mass regime of globular clusters (GCs, a star cluster in a tidal field reduces the M/L ratios of star 
Af < 2 X 10^ M0), dynamical M/L ratios tend clusters. There, we constrained the ranges of dissolution 
to be some 25% below predictions from simple timescales necessary for this loss of low-mass stars to quan- 
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titatively account for the drop of M/L observed for GCs. 
In the case of the Galactic GC system, it was found that 
dissolution timescales in the range = 0.6 — >20 Myr (cor- 
responding to total disruption times of i^"*'^' ~ 3 — 100 Gyr 
for a 10^ Mq cluster) are required to explain the observed 
M/L ratio decline. It was also shown that the M/L ratio 
decrease is strongest for low-mass GCs, which explains the 
observed c orrelation of increa s ing M/L ratio wi t h mas s dis- 
covered bv lMandushev et all (|l99lt) . iKruiissenI (|2008t) con- 
cluded that the scatter around this relation is caused by 
spreads in metallicity an d dissolution timesca le. 

As noted already in iKruiissen fc LamersI ([2008, here- 
after KL08), the next step is to apply these analytical clus- 
ter models including preferential loss of low-mass stars to 
individual clusters. This would then account for variations 
in metallicity and dissolution timescale. Such a study will 
naturally be restricted to GCs with measured M/L ratios 
for which realistic estimates of their individual dissolution 
timescale are available from information on their actual or- 
bit within the Milky Way potential. With the database of 
individual dissolution time scales at hand, the loss of low- 
m ass stars can be q uantified according to the prescriptions 
of IKruiissenI (|2008t) and KL08, leading to predictions for 
the drop of M/L for individual GCs. Those predictions are 
to be contrasted with the actual observed M/L ratios of 
these GCs. This will allow us to quantitatively test the hy- 
pothesis that the loss of low-mass stars is responsible for 
the too low M/L ratios of GCs, and hence also partially 
for the discrepancy of M/L between GCs and UCDs. 

Previous studies assessing the preferential loss of low- 
mass stars in Galactic GCs focus both on observations (e.g., 
De Marchi et al.l [20071) and theor y fe.g.. iBaumgardt et al.l 



20081) . In lDe Marchi et al.l (|2007D . the slopes of the stellar 



for dy- 
that of 
by the 



mass functions in GCs are measured for stars between 0.3 
and 0.8 Mq, thereby directl y reflecting possib l e low- mass 
star depletion. The study bv IBaumgardt et all (j2008( ) pre- 
dicts the same slopes using A''-body models and different 
degrees of mass segregation, assuming dissolution by two- 
body relaxation. The aforementioned papers both do not 
consider the M/L ratios of the GCs in question. 

In this study, the reference sample 
namical M/L ratios of G alactic GCs is 
iMcLaughhn k van der Marell (j2005f ). obtained 
fitting of single-mass King profiles. It contains data for 
38 GCs. Only a subsample can be used for our analysis, 
namely those clusters for which accurate proper motions 
and radial velocities are measured and can be translated 
to orbital parameters. Next to the destruction rates due 
to two-body relaxation, also the influence of disc shocking 
on the cluster dissolution needs to be taken into account. 
Both have to be derived from the orbital parameters. 
Several studies in which orbital information is derived 
and used to co mpute destruction r ates ar e available in 
the literature ( Gnedin fc Ostrikeil 119971 : iDinescu et al.l 
[l99llAllen et alll2006Ll2008ir all with certain beneflts and 
trade-offs. Specifically, iGnedi n fc Ost rike^ (|1997D assign 
statistically sampled orbits conforming to the bulk motion 
of the GC system in an axisymmetric pote ntial to derive 
destru ction rates of 119 globular clusters. IDinescu et al.l 
(|1999| ) use proper motions and radial velocities to compute 
the orbits and destruction rates o f 38 clusters in two 
axisymmetric potentials (Paczy nskil 119901 : iJohnston et al.l 
ll995f) . lAllen et all (|2006l . 120081) follow the same procedure, 
but consider both axisymmetric and barred potentials 



(jAllen fc Santillanlll991l:lPichardo et al.ll2004 respectively) 
for 54 globular clusters. They do not find a significant 
deviation between the results for both potentials. However, 
they do note that their calculated destruction rates are 
multiple orders of magnitude smaller than others in litera- 
ture and recommend combining their orbital information 
wi th the more ri gorous Fokker-Planck approach used by 
I Gnedin fc Ostriker (1997) to derive the destruction rates. 

We choose to adopt the or bital data and destr uction 
rate due to disc shocking from IDinescu et all (jl999[ ). Our 
study cannot be based on statistically assigned orbits but 
requires the actual orbits of in dividual clusters, thus exc lud- 
ing the estimated orbits from, Gnedin fc Ostrikej (fl99l . In 
addition, the IDinescu et al.l (|1999D destruction rates seem 
to be in better agreement with the observations than those 



from lAllen et al.l (2006, 2008). 

In Table [T] the observed properties are listed of the 24 
Galactic globular clusters for which the T^-band mass-to- 
light ratios (M/Ly) and orbital par ameters are available , 
i.e., the sample that is covered both by Dinescu et al.l (|l999t ) 
and lMcLaughlin fc van der Marell (|2005D . A first inspection 
of the observed M/Ly ratios can be made by comparing 
them to the ('canonical') M/L ratios from SSP models, 
which only depend on metallicity due to the invariance of 
the shape of the stellar mass function in these models. In 
Fig. [T] (left), the observed M/L ratios of our sample are 
plotted versus the canonically expected values that were 
computed by interpolating SSP models. These were emu- 
lated with the models from KL08 neglecting the preferen- 
tial loss of low-mass stars. The discrepancy between ob- 
served and expected M/Ly ratio is evident, as the canon- 
ical M/Ly are constrained to a much narrower and gen- 
erally higher range than the observed ones. The number 
histogram of the two M/Ly (Fig.[Tl right) further substan- 
tiates this dissimilarity. The observed M/Ly ratios are on 
average 74ly% of the canonically expected values. 

With the present paper we aim to quantify the con- 
tribution of dynamical effects such as the preferential loss 
of low-mass stars and the selective loss of stellar remnants 
(see KL08) to the discrepancy between the observed and 
canonically expected M/L ratios. In Scct.|2l we summarise 
the cluster models from KL08 and highlight the aspects 
that are particularly relevant to this study. The dissolution 
timescales for the cluster sample are computed in Sect. |3l 
whereas the predicted mass and M/Ly evolution are con- 
sidered and compared to the observations in Sect. H] We 
predict slopes of the low-mass stellar mass function and 
discuss observational tests to verify the preferential loss 
hypothesis for appropriate clusters in Sect. [5] In the final 
Sect. |6l we discuss the results and present our conclusions. 



2. Cluster evolution models and M/Ly evolution 

In order to study the evolution of clusters on specific or- 
bits, we use analytical cluster models (SPACE, see KL08) 
that incorporate the effects of stellar evolution, stellar rem- 
nant production, cluster dissolution and energy equiparti- 
tion. They are summarised here and are treated in more 
detail in KL08. In the second part of this section, the de- 
pendence of mass-to-light ratio evolution on initial mass, 
metallicity and dissoluti on timescale is a ssessed (for a more 
detailed description, see lKruiissenll2008l ) . 
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Cluster properties 
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Table 1. Observed properties for the cluster sample. Consecutive columns list the cluster NGC number, logarithm of the 
present-day cluster mass M, observed y-band mass-to-light ratio {M/ Lv)ohs, metallicity [Fe/H], galactocentric radius 
Rgc, King parameter Wq and canonically expected V^-band mass-to-light ratio (Af/Ly)can- Masses and observed M/Ly 
ratios are dynamical values. For all clusters, the standard error in [Fe/H] is assumed to b e 0.15 when computing th e 
error propagation (see Sect.[3|), which represents a conservative accuracy estimate (see, e.g.. ICarretta fc Gratton|[l997[ ). 
These errors determine the uncertainty on (A//iy)can in the last column, since the canonical M/L ratio only depends 
on me tallicity. 
*Prom 
I^From 



McLaughlin fc van der Marell (|2005D . 

Harris! ( 19961 ). The Rrrr. values used to compute the orbits in lDinescu et al.l ()1999D are from Zinn, private commu- 



nication. In extreme cases this m ay cause a sm all disagreement between the galactocentric radius quoted here and the 
apogalactic distance predicted bv lDinescu et al. (1999) (see Table 
*The value for NGC 5139 (cjCen) is derived from lBedin et al.l (|200l . 



2.1. Summary of the models 

In the SPACE cluster models, clusters gradually lose mass 
due to stellar evolution and dissolution. The total cluster 
mass evolution is determined by 



dAfci 
dt 



fdM, 



\ dt 



fdM, 



\ dt 



(1) 



dis 



where the first term denotes mass loss due to stellar evolu- 
tion and the second represents mass loss by dissolution. 
Additionally, the formation of stellar remnants and the 
mass-dependent loss of stars by dissolution are taken into 
account, thus providing a description of the changing mass 
function and cluster mass in remnants. 

Stellar evolution is included by using the Padova 
1999 isochronesQ. These are available for metallici- 
ties Z = {0.0004,0.004,0.008,0.02,0.05} (or [Fe/H] = 
{-1.7, -0.7, -0.4, 0.0, 0.4}), which thus restricts our model 
computations to these values. Stellar evolution removes 
the most massive stars from the cluster and increases the 



^ These isochrones are based onlBertelli et al.l (|1994| ). but use 
the AGB treatment as in iGirardi et al.1 (|2000t ). 



non-luminous cluster mass by turning stars into remnants, 
which is included by adopting an initial-remnant mass re- 
latioifl. A Kroupa (2001) IMF is assumed. 

Cluster dissolution represents the dynamical cluster 
mass loss due to stars passing the tidal radius, which acts 
on the timescale Tdis: 



/ dAfei(t) 
I dt 



M,i{t) _ M,i{t) 



1-7 



dis 



Tdi: 



to 



(2) 



where Mc\{t) represents the present day cluster mass and 
the second e quality follows from the relation derived by 
iLamers et al.i (2005a. ): 



Tdii 



to[M,iit)/MQ] 



(3) 



The characteristic timescale to depends on the environment 
and determines the strength of dissolution. For example. 



For white dwarfs, this relation is taken fro m iKalirai et al 



( 200^ ), while for neutron stars the relation from lNomoto et al 



^9881) is used. Black hole masses are a ssumed to be c onstant at 
8 Mo , in agreement with observations (|Casaresll2007D . For more 
details, see KL08. 
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Fig. 1. Left: Observed mass-to-light ratio [M / Lv)ohs versus the canonically expected mass-to-light ratio {M / Lv)ca,m 
together with their la standard errors (see Table [T]). The dotted line follows the 1:1 relation. Clusters for which the 
disagreement is larger than la are plotted as dots. Right: Number histogram of {M/Lv)can (dia monds, shaded area) and 
(M/ Ly)obs (tr iangles, hashed area). For comparison, canonically expected M/Ly ratios from iBruzual fc CharlotI (|2003t ) 
for a lChabrierl ((2003) IMF are overplotted (squares, dotted hue). Again, the error bars denote la deviations, which were 
determined from 30k random realisations of the underlying data. 



in the case of dissolution by two-body relaxation, to de- 
pends on tidal field strength and therefore on the angular 
velocity of t he cluster orb it. Typical values are to = 10^ — 
10^ yr (e.g.. iLamers et al . 2005b), translating into a total 
disruption time t^^l'^^ « 10^—10" yr for a 10^ Mq cluster. 
The expon ent 7 is found to be 7 w 0.62 , both from ob - 
servations (iBoutloukos fc Lamerg|[2 003: Gi eles et"al]|2005[ ) 
and from the 'Baum gardt fc Makind (|20030 iV-body simu- 
lations of tid al dissolution for clu sters with King parame- 
ter Wo = 5 fLamers et al . 2005b). However, it is recently 
derived by iLamers et al.l (2009,) that 7 = 0.70 for King 
parameter Wq = 7. Since this concentration more closely 
resembles the mean King parameter for Galactic GCs (see 
Table[T]), we adopt 7 = 0.70 throughout this study. 

The effect of dissolution on the mass function depends 
on the dynamical state of the cluster. As it evolves towards 
energy equipartition, low-mass stars are preferentially lost 
from the cluster. This mass loss (in the 'preferential mode', 
KL08) is approximated b y increasing the minimum stellar 
mass ()Lamers et al.l [20061 ) . while evaporation that is inde- 
pendent of stellar mass (mass loss in the 'canonical mode', 
KL08) is accounted for by decreasing the normalisation of 
the mass function. In our models, both modes coexist to 
allow for intermediate modes of mass loss. Their relative 
contributions are fitted such that th e M/ Ly ratio evolution 
match es the A^-body simulations bv lBaumgardt fc Makind 
(|2003l) . 

Cluster photometry is computed by integrating the stel- 
lar mass function over the stellar isochrones, yielding clus- 
ter magnitude evolution M\{t, Mdj) for a passband A and 
a cluster with initial mass Md.i. 



2.2. Dependence of mass-to-light ratio on model parameters 

The models described in Sect. 12.11 vield a mass-to-light ra- 
tio evolution that depends on the dissolution timescale, 
metallicity and initial cluster mass. In canonical models, 
i.e., without the preferential loss of low-mass stars, M/L 
monotonously increases with time. For a given age and 
metallicity, these models provide M / L ratios that are in- 
dependent of cluster mass. On the other hand, our mod- 
els including dynamical effects predict a mass-dependent 
drop in mass-to-light ratio^due to the ejection of low-mass, 
high-M/L stars ()KruiissenH2008L KL08). In Fig. El the V- 
band mass-to-light ratio evolution M/Ly is shown for two 
metallicities and several initial cluster masses. In both pan- 
els, the upper curve marks the canonical mass-to-light ratio 
evolution, while the others correspond to cluster evolution 
including the preferential loss of low-mass stars for differ- 
ent initial masses. Since low-mass clusters evolve on shorter 
timescales than massive ones, the deviation of their mass- 
to-light ratio evolution with respect to canonical models 
arises at earlier times than for massive clusters. 

The mass-to-light ratio decrease can be quantified by 
considering the ratio of the observed or predicted M/ L to 
its canonical value to divide out their metallicity depen- 
dence: 

^ ^ (M/j^y)obs/prcd . . 

{M/ Lyjcan 

Figure [31 shows the predicted fraction of the canonical 
M/L ratio Qprcd as a function of i/i^"***' for clusters with 
initial masses in the range Mc\;i — 2 x 10^ — 10^ M©, 
dissolution timescales to = {1: 10} Myr and metallicities 
Z = {0.0004,0.004}. It shows that Qprcd is independent 
of metallicity, initial cluster mass and dissolution timescale 
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Fig. 2. Left: y-band mass-to-light ratio evolution M/Ly 
for to — 1 Myr, Z — 0.0004 and initial masses in the range 
logMci^i = 5 — 8 with G.25-dex intervals. Right same graph, 
but for Z — 0.004. From top to bottom, different curves 
represent the M/ Ly evolution for decreasing initial cluster 
masses. 



when considered as a function of the elapsed fraction of the 
total disruption time t/t^^^l'^^. The three-component linear 
approximation illustrates the well-defined uniform correla- 
tion and is given by 



1 



Qprcd — 



for t/t^°l''^ < 0.2, 
1.142 - 0.71i/t*d°f' for 0.2 < t/t*d°s''' < 0-7, 
1.471 - 1.18t/t*d°f' for t/t^Jll''^ > 0.7, 



(5) 



which applies for all initial conditions, i.e., is independent 
of the cluster properties or environment Equation [5] does 
not include possible effects of primordial mass segregation 
on the change of the mass function. From model runs where 
we did assume the preferential depletion of low-mass stars 
from t = on we know that its effects become about 10% 
stronger with respect to purely dynamically induced low- 
mass star depletion (KL08). This number should be treated 
with some care, because our models are based on A^-body 
simulations of cl usters that did not start ou t in a mass- 
segregated state (jBaumgardt fc Makinoll20d3[ ). 

The relation between the predicted fraction of the 
canonical mass-to-light ratio Qprcd and the elapsed fraction 
of the total disruption time is expected, since in our models 
the decrease of M/L is the result of dynamica l evolu tion. 
It is in agreement with studi es by [Richer et al.l (|l99lD and 
iBaumgardt fc Makinol (|2003( ). who find that the depletion 
of the low-mass stellar mass function in globular clusters 
is closely related to the elapsed fraction of the total dis- 
ruption time. Considering the physical processes driving 
dissolution, the result is not surprising either. Two-body 
relaxation is known to preferentially eject low-mass stars 
(jHenon|[l969t ) and tidal shocks remove the outer parts of 
the cluster, which in the case of mass segregation are con- 
stituted by low-mass stars. 

The evolution of cluster mass and M/Ly can both be 
considered in the {M, M/Ly}-plane. The resulting 'evolu- 
tionary tracks' are shown in Fig. 3] for two different dissolu- 
tion timescales and again for two metallicities and a range 
of initial cluster masses as in Fig. [51 Clusters start with their 




1.2 



Fig. 3. The ratio of predicted to canonical mass-to-light 
ratio Qprcd as a function of the elapsed fraction of the total 
disruption time i/idJ**^'. Dotted curves denote model predic- 
tions for a broad range of initial conditions (varying initial 
masses, dissolution timescales and metallicities), while the 
solid line describes a three-component linear approximation 
to the models. 



4 - 



Z = 0.0004 
!„= 1 Myr 

• t = t„ 



pref 




Z = 0.0004 
to=10 Myr 



^ Please note that a Kroupa IMF was assumed here. For sub- 
stantially different IMFs the relation will vary. 



7 8 5 6 
log Mass [Mg] 

Fig. 4. Cluster evolution in the {Af, M/Lyj-plane for 
to = {1, 10} Myr and Z = {0.0004, 0.004}. Solid curves rep- 
resent cluster evolutionary tracks for initial cluster masses 
in the range Mc\;i = 10^ — 10^ M0with 0.5-dex intervals. 
Cluster isochrones at i = 12 Gyr are described by dashed 
lines, while these at t = 10 and i = 14 Gyr are denoted by 
dotted fines (bottom and top, respectively). Dots denote 
the onset of the preferential loss of low-mass stars for each 
evolutionary track. 



initial masses and with M/Ly ratios close to zero, corre- 
sponding to an initial position on the x-axis of Fig. |4l As 
time progresses, clusters initially evolve to lower masses and 
increasing M/Ly due to the death of massive stars, trans- 
lating into up- and leftward motion in the {M^M/Ly}- 
plane. When the preferential loss of low-mass stars becomes 
an important mechanism (the onset of which is marked by 
dots for each evolutionary track), the M/Ly increase is 
turned into a decrease instead, as also illustrated in Fig. O 
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Orbital parameters 
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Table 2. Orbital parameters for the cluster sample, together with their Icr standard errors. Consecutive columns list 
the cluster NGC number, apogalactic radius i?a, perigalactic radius Rp, eccentricity e, orbital period P, circular velocity 
of the gravitational potential at the distance of apogalacticon Vc^a. and destruction rate due to disc shocking Vgh- The 
circular velocities are computed in the galactic plane (z = 0). Because the gravitational potentials of the disc and bulge 
decrea se with \z\, thi s impli es that for clusters with < 10 kpc the actual Vc^a can be 5 — 15% lower. 
*From lDinescu et al.l (|1999 D. 

^Computed using the galactic potential from lPaczvnskil (|l99Cll ). 



In Fig. m the thus attained maximum in the M/Ly evolu- 
tion is best visible for low cluster masses and to = 1 Myr. 

Since Galactic globular clusters generally share the same 
ages (e.g.. IVandenberg et al.llT990l ). the observed distribu- 
tion of GCs in the {M, M/Ly}-plane would follow curves 
of equal age in Fig. [5] if there were no spreads in metal- 
licity and dissolution timescale. These curves, or cluster 
isochrones, are shown for ages t = {10, 12, 14} Gyr. Along 
the isochrones, M/Ly increases with cluster mass since 
massive clusters have spent a smaller fraction of their total 
disruption time than low-mass clusters and will therefore 
have experienced a smaller M / Ly decrease due to low-mass 
star depletion. The curves flatten at the highest masses, 
since these clusters have not yet exhibited significant pref- 
erential low-mass star ejection. 

From Fig. [5] we infer the influences of dissolution 
timescale and metallicity on the mass-to-light ratio evo- 
lution. The dissolution timescale sets the cluster mass for 
which the down-bend of the cluster evolutionary tracks can 
occur and therefore also determines the location of the 
'knee' in the cluster isochrones. The metallicity determines 
the vertical extent of the cluster evolutionary tracks and 
thus the M /Ly-normalisatio n of the cluster isochrones. As 
set forth in iKruiissenI ()2008f ). the natural spread in disso- 
lution timescale and metallicity thus explains the scatter 
around t he relation between M / L and cluster mass ob- 
served by iMandushev et "all ()1991[ ). 



3. Determining the dissolution timescale 

To assess the influence of the preferential loss of low-mass 
stars on the low observed mass-to-light ratios, the orbital 
parameters of individual clusters are to be translated into 
the appropriate dissolution timescales to for use in our clus- 
ter models. The computation is treated in this section. 

3.1. Cluster dissolution timescales from orbital parameters 

For globular clusters, dissolution due to two-body relax- 
ation in the Galactic tidal fielcfl an d disc shocking are th e 
main dissolution mechanisms (e.g.. [Chernoff et al.lll986l ). 
The total dissolution timescale to, tot can be written as 



^O.tot ^O.cvap ^O.sh 

where to,cvap denotes the dissolution timescale due to two- 
body relaxation or evaporation (carrying the subscript 
'evap') and t^ sh the dissolution timescale due to disc shock- 
ing. 

For the dissolution timescale due to two-body relax- 
ation, we use the expression for the total disruption time 



* This includes the effect of bulge shocks, which occur in clus- 
ters on eccentric orbits. 
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fro mlBaumgardt fc Makinol (|2003l . Eq. 10) as approximated 
bv lLamers et all (|2005bf) to write 



Dissolution timescales (7 = 0.70) 



^220 km s 



-1 



-1 



(1-e), (7) 



with ijfcvap dissolution timescale due to two-body re- 
laxation for a circular orbit at the solar galactocentric 
radius, i?gc,a the apogalactic radius of the cluster orbit, 
Vca the circular velocity of the gravitational potential at 
the distance of apogalacticon and e the orbital ec c entric - 
ity. Values for i?gc,a are taken from iDinescu et aT] (|l999t) . 
while the circular velocities a r e com puted for the galac- 
tic potential from iPaczvnskil (|1990D. This potential, as 
well as the one fro m iJohnston et al.l (|1995D . is used by 
iDinescu et al.l (|l999l ) in the determination of the cluster or- 
bit s. By comparing our models to the A^-body simulations 
bv lBaumgardt fc Makinol (|2003[ ) we find tf ,,^^p = 21.3 Myr 
for clusters with Wq = 5 King profiles, in ve ry close agree- 
ment w ith earlier report ed values of 20.9 Myr ((Lamers et al.l 
l2005a^ and 22.8 Myr (jLamers fc GielesI I2006D . Using the 
same method for 7 = 0.7, corresponding to Wq = 7 King 
profiles (see Sect. [23]), we obtain <§^evap = lO-^ Myr. This 
is the adopted value in this paper. 

The dissolution timescale due to disc shocking can be 
obtained from the globul ar cluster dest r uction rates due 
to disc shocking Vgh from IDinescu et all (|l999t ). Following 
from Eq. [31 a present destruction rate vit) is related to a 
dissolution timescale by 



10 
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to(A/ciW/M0)7^ 



(8) 
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Table 3. The dissolution timescales (for 7 = 0.70) due 
to two-body relaxation (io,cvap), disc shocking (io,sh) and 
both mechanisms (io,tot), together with their la standard 
errors. All values are rounded to one decimal. 



with ly in units of (10 Gyr)^^, to in years, and Mci{t) de- 
noting the cluster mass at age t. The denominator repre- 
sents an estimate for the total cluster lifetime. This expres- 
sion can be inverted to obtain to.sh from Vsh- However, in 
IDinescu et al.l (|l999t ) constant M/Ly = 3 M0 Lq^ is as- 
sumed to compute the cluster masses. Since their destruc- 
tion rates are derived from a relation i^sh oc M^^, these 
should be corrected for the actual mass-to- light ratios. We 
define the correction factor 



{M/L 



V)cst 



{M/L 



V)ohs 



(9) 



constant mass-to-light ra- 
10 Lq^ and the denomina- 



with the numerator the 
tio {M/Lv)cst = 3 Mg 
tor the observed dynamical ma ss-to-light ratio from 
[McLaughlin fc van der Marell (|2005f) (see Table [1]). This al- 
lows us to express the dissolution timescale due to disc 
shocking as 



10 



10 



iQ.sh 



■^corr^sh 

(Mei(t)/M0)'r- 



(10) 



Substitution of Eqs. [7] and [TO] into Eq. [5] then yields the 
total dissolution timescale tg.tot- 

In Table [2] ou r cluster samp l e is lis ted with the orbital 
param eters from IDinescu" et al.l (|l999t ) for the IPaczvnskil 
(|199 0^ potential and our computed circular velocities of 
the gravitational potential at the distance of apogalacti- 
con. The corresponding dissolution timescales can be found 
in Table [3] The values for the dissolution timescale range 



from io.tot = 0.5 — 20 Myr, corresponding to total disrup- 
tion times for a 10^ M0 cluster in the range t^'?*'*' = 8 — 
300 Gyr. This is is in good agreement with the range that 
is required for low-mass star de pletion to explain the ob- 
served mass-to-light ratio drop (jKruiissenl [20081. By com- 
paring the dissolution timescales for two-body relaxation 
^o.cvap and disc shocking to.sh, we can see that the latter 
destruction mechanism is important (i.e., lowers the total 
dissolution timescale io.tot by more than 40% with respect 
to to,evap) for the clusters NGC 288, 5139 (wCen), 6362, 
and 6779. These clusters all have perigalactic radii smaller 
than 3 kpc (see Table [5]) . For the error analysis of Tables [5] 
and [3] and of the rest of this paper we refer to Appendix [X] 



4. Predicted and observed mass-to-light ratios 

In this section we combine our cluster models and the de- 
rived dissolution timescales to study the mass-to-light ra- 
tio evolution for our sample of 24 Galactic globular clusters. 
Present-day M/Ly ratios are predicted for the cluster sam- 
ple and are compared to the observations. We also discuss 
the possible causes for the individual clusters that still lack 
convincing agreement. 



4.1. Predicted mass-to- light ratios for the cluster sample 

We employ the cluster models treated in Sects. [2] and [3] to 
predict M/Ly ratios for the cluster sample. The input pa- 
rameters for the models are the dissolution timescale to and 



8 



J. M. D. Kruijssen and S. Mieske: On the reduced mass-to-light ratios of Galactic globular clusters 




Fig. 5. Top left: Time-evolution of M/Ly for NGC 5466. 
The solid line represents the Z — 0.0004 model with the 
dissolution timescale of the cluster {to = 10.2 Myr), while 
the dotted curve indicates the canonical M/Ly evolution, 
i.e., if the preferential loss of low-mass stars were omitted. 
The dashed line denotes constant age of t = 12 Gyr. The 
predicted M/Ly of NGC 5466 is marked by a cross and the 
onset of the preferential loss of low-mass stars is specified 
with a dot. Top right: Evolution in the {M, M/Ly}-plane 
for NGC 5466. Curves and symbols have the same meaning 
as in the top-left panel. Bottom left: same as top left, but 
for NGC 6779 (to = 1.2 Myr). Bottom right: same as top 
right, but for NGC 6779. 



4.0 4.5 5.0 5.5 6.0 6.5 7.0 
log Mass [Mq] 

Fig. 6. Observed (crosses) and predicted (triangles) dis- 
tribution of GCs in the {A/, Af/Ly}-plane. The solid line 
represents our linear fit to the observati ons, while the dot- 
ted li ne denotes the relation found by iMandushev et all 
([il9il).The error bars in the top left corner denote the av- 
erage la uncertainty on the observations, while the error 
bars in the bottom right corner represent the average la 
uncertainty on the predictions. 



metallicity Z. The latter is derived from the iron abundance 
[Fe/H] (see Table [J) according to 

Z = Zq X 10P°/«1, (11) 

with Zq = 0.02, while the dissolution timescale is taken 
from Table [3] (to,tot)- Since all clusters in the sample have 
metallicities Z < 0.004, for each cluster the models are 
computed with metalhcities Z = {0.0004,0.004} and the 
appropriate dissolution timescales. The evolution is com- 
puted for a grid of initial cluster masses, yielding cluster 
evolution tracks for the mass and y-band mass-to-light ra- 
tio M/Lv- For both metallicities, at i = 12 Gyr the tracks 
are interpolated over the mass grid to match the observed 
cluster mass. This provides predictions for M/Ly, the ini- 
tial cluster mass Mdj and the total disruption time t^"*^' for 
two metallicities. These are then interpolated over metal- 
licity to obtain the model predictions for the appropriate 
metallicity. 

Examples of the M/Ly evolution with time and mass 
are shown in Fig. Ofor NGC 5466 and 6779. The predicted 
M/Ly are slightly offset with respect to the model curves 
because the models here are computed at Z = 0.0004 while 
the predictions are interpolated over metallicity. However, 
the variation with metallicity is small for the displayed 
clusters, since their metallicities are close to Z = 0.0004. 
It is evident that low-mass star depletion has a much 
stronger effect in the case of NGC 6779 than for NGC 5466. 
Considering their dissolution timescales (to = 1-2 Myr ver- 
sus to = 10.2 Myr, respectively) and the resulting mass 
evolution, this is not surprising since NGC 6779 has suf- 
fered much stronger mass loss than NGC 5466. 



The predicted mass-to- light ratios M/Ly, initial masses 
Mcij and remaining lifetimes tj^°*''' — t are listed for our en- 
tire GC sample in Table [H In addition, the observed and 
predicted fractions of the canonical M/Ly ratios Qobs and 
Sprcd are shown, as well as the agreement between our pre- 
dicted M/Ly and the observed values. Combining Tables [3] 
and m we see that GCs with short dissolution timescales 
indeed have low predicted M/L ratios. 

4.2. Comparison of predictions to observations 

The fifth column in Table [3] indicates the ratio be- 
tween observed and predicted mass-to-light ratio Qobs = 
{M/Ly)ohs/{M/Ly)can- Analogously, the sixth column 
gives the ratio between predicted and canonical mass-to- 
light ratio Qprcd = {M/ Ly)p^ed/{M/ Ly)ca.n- On aver- 
age, the former ratio is 0.74tQ Q7, while the latter ratio is 
0.85±0.01 for the 24 GCs investigated. There are factors 
that introduce biases when comparing the predictions to 
the observations. Specifically, the observations are likely bi- 
ased to central M/L ratios for some GCs, while we predict 
global values. In Sect. 14.3] a more detailed consideration is 
provided in which the comparison of the predictions to the 
observations is refined. 

The seventh column in Table S] gives the level of agree- 
ment between the observed and predicted mass-to-light ra- 
tios, which is defined as n if (n — l)a < \AM/Ly\ < na for 
AM/Ly = {M/Ly)prcd-{M/Ly)ohs- Within the la uncer- 
tainty, the predicted M/ Ly agree with the observed values 
for 12 clusters out of the 24-cluster sample. A Gaussian 
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Model predictions 
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Table 4. Predicted T^-band mass-to-light ratio {M/ Lv)pred, logarithm of the initial cluster mass M^^i and remaining 
lifetime t^^^l^^ — t for the clusters under consideration, together with their la standard errors. In the fifth and sixth column, 
respectively the ratios of observed to canonical mass-to-light ratio Qobs and predicted to canonical mass-to-light ratio 
Qpred are listed. The seventh column gives the level of agreement bet-ween the observed mass-to-light ratio {M/Lv)ohs 
and predicted mass-to-light ratio {M/Lv)pT:cd, which is defined as follows. For AM/Ly = [M / Lv)i,icd ~ [M / Lv)ohs, an 
agreement value of "1" means |AM/Ly| < Icr, "2" means Icr < |AM/Ly| < 2tT, etc., with cr^ = -|- ctJm/l) 



distribution of errors would yield an expected 16 out of 24 
clusters to be found within Icr. 

As a first comparison and a nalo gously to the presenta - 
tion in i Mandushev et al.l (|l99lh and lReikuba et al.l (|2007f) . 
in Fig. [S] the distribution of GCs in the {M, M/Ly}-plane 
is shown for the observed and predicted mass-to-light ra- 
tios. Both populations fall within the same range and follow 
comparable trends of in creasing MjLy with cluster mass. 
IMandushev et al.1 ()1991[ ) already provided an expression for 
the observed logarithm of the mass as a function of mag- 
nitude, which allows for a derivation of the expected trend 
in Fig. m They fit 

log (M/Mq) = (-0.456 ± 0.024)My -f (1.64 ± 0.21), (12) 

where My represents the y-band absolute magnitude of 
the cluster. Adopting a solar value of My.c--) = 4.83, the 
relatio n between MjLy and mass from IMandushev et all 
(|1991l) can then be expressed as 

\o%{MlLy) = (-O.12±O.O5)log(M/M0)-(O.49±O.21).(13) 

A first-order Taylor expansion of MjLy around 
log (M/Mq) = 5.2 then gives 

MjLy « 0.38 log (M/Mq) - 0.55, (14) 

which has a linear slope of 0.38. The best fitting slope for 
our sample is 0.41 ± 0.28 , thus agreeing with the value from 
IMandushev et all (|l99ll ). The large uncertainty arises from 
the scatter in Fig. [HI 



The trend of increasing mass-to-light ratio with mass 
is expected from the models shown in Fig. |¥| and the dis- 
cussion in Sec. 12.21 However, there the slope is ^ 0.6 — 1.0 
for metallicities Z = 0.0004 — 0.004 and increases with Z. 
For some metallicities, the model slope is thus more than 
Icr steeper than the fitted slope. This is not surprising, be- 
cause the models each have a single dissolution timescale 
and metallicity, while in reality both quantities have a 
spread that causes horizontal and vertical scatter, respec- 
tively. It turns out that the spread in dissolution timescale 
has a stronger effect on M /L than the spread in metallicity 
(Kr uiissen.2008 ), implying that the scatter in the horizon- 
tal direction is largest and that the slope fitted to the entire 
sample is shallower than that of a single model. 

In Fig. |7l a more specific comparison is made between 
the observations and model predictions using the same 
framework as for the canonical expectations in Fig. [TJ 
Again, the left-hand panel plots the observed versus the 
predicted mass-to-light ratios, while the right-hand panel 
shows the number histograms of the two. In the left-hand 
panel it is shown that the predictions for half of the clusters 
are such that they reach down to the appropriate mass-to- 
light ratios. When comparing this panel to its analog in 
Fig. [Tl the improved agreement with the observations is ev- 
ident. Nonetheless, there is an aggregate of deviating GCs 
below the 1:1 relation at {M / Ly)p^cd ~ 2 M© Lq^, rep- 
resenting the clusters for which no strong low-mass star 
depletion is expected from the models due to their long 
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Fig. 7. Left Observed mass-to-light ratio {M / Lv)ohs versus the predicted mass-to-light ratio {M/ Lv)prcd, together with 
their la standard errors. The dotted line follows the 1:1 relation. Clusters for which the disagreement is larger than la 
are plotted as dots. Right: Number histogram of {M / Lv)picd (diamonds, shaded area) and {M/Lv)ohs (triangles, hashed 
area). Again, the error bars denote la deviations, which were determined from 30k random realisations of the underlying 
data. 



disruption times. Consequently, the predicted M/Ly for 
these clusters are similar or equal to the canonical values. 
Except for NGC 6809, there are no clusters above the 1:1 
relation that are inconsistent with the observations. The 
number histogram of the observed and predicted mass-to- 
light ratios in the right-hand panel of Fig. [7] confirms both 
the improved agreement between observed and predicted 
M/Ly with respect to Fig. [1] and the accumulation of a 
number of clusters near the canonical M/Ly in the model 
predictions. 

4.3. Discussion of discrepant clusters 

In total, there are twelve clusters with a worse than la 
agreement between the model predictions and observations. 
Five of these have worse than 2a agreement, while we would 
expect only one. Here, we discuss possible reasons behind 
the discrepancy. 

The deviant clusters below the 1:1 relation in the left- 
hand panel of Fig. [1 being NGC 104, 1851, 1904, 4147, 
4590, 5272, 5904, 6205, 6341, 6934 and 7089, generally share 
properties such as relatively wide orbits and long dissolu- 
tion timescales. Due to their long dissolution timescales, 
they are all predicted to have near-canonical M/Ly. This 
is illustrated in Fig. [51 where the fraction of (M/iy)obs 
and (M/Ly)pi.od with respect to the canonical {M / Ly)ca,n 
is shown in panels similar to Fig. [T] In both the left- and 
right-hand panels of Fig. [51 the accumulation of too high 
predicted mass-to-light ratios occurs near or at Qpred — 1- 
Since per definition (M/Ly )prcd < (Af/Ly)can, no values 
Qprcd > 1 are found. In that range, the apparent disagree- 
ment between the observed and predicted histograms is dis- 
putable since all but one cluster (NGC 6809) are in la- 
agreement with their canonical mass-to-light ratios. 



While our predicted mass-to-light ratios are 
global (i.e., cluster- wide) values , the observations 
from [McLaughlin fc vaii der Marell (|2005[) are de- 
rived from central v el ocity dispersion measurements 
from iPrvor fc MevlanI (|1993f ) and are extrapolated 
to global values using s ur face brightness profiles 
(McLaughli n fc van der Marell l2005l and references 
therein). They fit isotropic single-mass King models and 
thus neglect any radial gradients of M/L ratio or mass 
function slope. Consequently, the values of (M/I/y)obs do 
not contain any information about such gradients and for 
some clusters only accurately reflect the M/Ly ratio in 
their central parts. The global M/Ly ratios of clusters 
with strong radial M/L gradients are at best approximated 
(McLaughlin, private communication). For instance, the 
centre of a mass-segregated cluster may be populated with 
massive, i.e., luminous stars, yielding a lower M/Ly than 
its global value. 

The disagreement between the global and central M/L 
is expected to be largest for clusters that have suffered 
relatively weak mass loss but are internally evolved. In 
that case, the low-mass stars are outside the core but 
still bound to the cluster and are included in the global 
M/L, w hile they do not play a role in th e value de- 
rived bv [McLaughlin fc van der Marell (|2005[) . This is in- 
deed the case for the discrepant GCs in our sample, which 
not only have long dissolution timescales but also higher 
King parameters Wq, implying that mass segregation can 
be reached on relatively shorter timescales. For the GCs 
with worse than la agreement below the 1:1 line in Figs. [7| 
and[8]we have average King parameter Wq — 7.6, while for 
the Icr-consistent GCs we find Wq = 6.1, both with stan- 
dard errors < 0.1. This further validates our explanation for 
the difference between the central and global M/L ratios 
of these GCs. 
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Fig. 8. Left: The ratio of the observed mass-to-light ratio {M/ Lv)ohs to the canonically expected {M/ Lv)can versus the 
ratio of the predicted mass-to-light ratio (A//L\/)prod to [M / Lv)cam together with their la standard errors. The dotted 
line follows the 1:1 relation. Clusters for which the disagreement is larger than Icr are plotted as dots. Right: Number 
histogram of the ratio Qprcd (diamonds, shaded area) and Qobs (triangles, hashed area). Again, the error bars denote Icr 
deviations, which were determined from 30k random realisations of the underlying data. 



In a recent study, iDe Marchi et al.l (j2007f ) find that ex- 
tended (low-concentration and Iow-Wq) GCs are depleted in 
low-mass stars, which they confirm to be in accordance with 
predictions by theoretical studies (IChernoff fc Weinberg! 
[1990; Takahashi & Portcgies Zwari|2OOO0, while GCs with 
high values of Wq have close to canonical mass functions. 
This is in agreement with our predicted M/L for these clus- 
ters and suggests that the observed M/Ly are indeed un- 
derestimated. A more precise check can be made b y com- 
paring the low-mass star depletion from Do Marc hi et al.l 
(|2007l ) with the observed and predicted fractions of the 
canonical M/Ly ratios Qobs/prcd- This can be done for four 
GCs with worse than Icr agreement, being NGC 104, 5272, 
6341 and 6809. For NGC 104 and 5272 the observed de- 
pletion is not strong enough to draw any definitive con- 
clusions, while f o r NG C 6341 and 6809 the results from 
iDe Marchi et al.l (|2007f ) are clearly more consistent with 
our predictions than with the observed M/Ly (see also 
Sect. [S] and Fig. [5]). This substantiates the claim that some 
GCs have observed M/L ratios that are biased to lower 
numbers. To test this assertion, global observational mea- 
surements of the velocity dispersion would be needed to 
enhance the accuracy of the present observed M/L ratios. 

We now revisit the mean M/Ly fractions of the canoni- 
cal value presented in Sect. l4.2l bv leaving out the GCs that 
may have strongl y different central and g l obal M /L ratios. 
It was shown by iBaumgardt fc Makinol (|2003f) that for a 
10^ Mq cluster core collapse is reached within a Hubble 
time if Wq > 7. This timescale is increased by a factor three 
for a GC with typical initial mass of 10^ Mq, but mass seg- 
regation manifests itself on a shorter timescale than the 
core collapse time. The relative mass loss due to dissolu- 
tion of a 10^ Mq GC is smaller than 10% after 12 Gyr for 
dissolution timescales fo.tot > 5 Myr. These limits could 
separate GCs with similar global and central M/L ratios 



from those with pronounced differences between the two. 
We exclude GCs with both to, tot > 5 Myr and Wq > 7, as 
well as NGC 6809 (which has a M/Ly ratio that cannot 
be explained by any model as it is 1.6 times the canoni- 
cal value). This yields an average observed fraction of the 
canonical M/Ly of Qobs — ^■'^^^^'ai ^'^d ^ predicted value 
of Qpred = 0.78 ± 0.02. For the excluded GCs, we have 
QXs' = 0.68t°:g^ and Q^^f^ = 0.96 ± 0.01, reflecting the 
fundamental difference between both values. Although the 
cuts we made represent only 'educated guesses', it is evi- 
dent that the agreement between theory and observations 
is much better for those GCs for which we can be more cer- 
tain that the central M/Ly reflects the global value. For 
these GCs, our models confirm an average M/L ratio drop 
of about 20% due to low-mass star depletion, corresponding 
to about 1/4 of the observed difference in M/Ly between 
GCs and UCDs. 

Another option could be that the dissolution timescales 
of GCs on wide orbi ts are overestimated (as suggeste d 
for different reasons bv lKruiissen fc Portegies Zwartll2009[ ). 
possibly due to a dissolution mechanism tha t is not included 
in our analysis. White dwarf kicks (Freg eau et al.l 120091 ) 
could be a candidate for such a mechanism. This would im- 
ply that some of our predicted dissolution timescales and 
M/L ratios are overestimated. Also, we do not assume clus- 
ters to be initially mass-segregated. Some of the clusters un- 
der consideration here are likely not to have reached energy 
equipartition within a Hubble ti me, but still exhibit ev i- 
dence for mass segregation (e.g., lAnderson fc Kinj|1996[ ). 
This points to prim ordial mass segregation in these cases, 
which is shown by IBaumgardt et alT (|2008f ) to effect ad- 
ditional low-mass star depletion that we did not account 
for (see also Sect. 12.21) . The additional modeled M/L ra- 
tio decrease would be ^ 10% (KL08). However, this is not 
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sufRcient to lift the discrepancy for any of the deviating 
GCs. 



5. Observational verification 

If the decrease of M/L ratio with respect to the canonical 
value is indeed due to low-mass star depletion, one would 
expect a correlation between the observed slope of the low- 
mass MF ttobs and the ratio of the predicted and canonical 
M/Lv ratios Qprcd- Specifically, for a powerlaw MF with 
n oc TO~", a small value of Qprod would be signified by a 

reduced value of Oobsj 

In a study by iDe Marchi et alj (|2007[ ). MF slopes are 
determined in the stellar mass range m = 0.3 — 0.8 M© for 
several Galactic globular clusters, based on a compilation 
of results fro m HST imaging of different sources. By re- 
analysing the iBaumgardt fc Makiiiol (|2003f ) TV-body data , 
iBaumgardt et all (|2008[ ) conclude that for a lKroupal (|200l[ ) 
IMF the canonical slope in that mass range is ao = 1.74, 
which is thus expected to be measured for clusters with 
canonical M/L ratios or Qpred = 1. In addition, they pro- 
vide a fou rth-order powerlaw fit to the iV-body simula- 
tions from I Baumgardt fc Makinol ()2003 f) for a as a func- 
tion of the elapsed fraction of the total disruption 
time. By inverting our relation between Qprcd and i/^dis''' 
(Eg. and inse r ting t he outcome into Q;(t/ij°*'*') from 
IBaumgardt et al.l (|2008l Eq. 4), we obtain an expression 
for the predicted MF slope aprcd between 0.3 and 0.8 M© 
as a function of the fraction of the canonical MjLy ratio 
Qprcd- Inversion of Eq. [5] yields 



total 



1.25 - 0.85Qpred for 
1.61-1.41Qprcd for 
0...0.2 for 



^prcd 



< 0.645, 



0.645 < Qprcd 

Qprcd = 1, 



<1, (15) 



where the uncertainty tjt^^^^^ = 0...0.2 for Qpred = 1 arises 
due to the range of t/^dis'^' over which it is constant in 
our models. Combinatio n of this expression and Eq. 4 from 
IBaumgardt et all (|2008f ) then provides the relation between 
apred and Qpred: 



^prcd 



^a„Q;jrcd for Qprcd < 0.645, 



n=0 
4 



I] ^" Qprcd for 0.645 < Qprcd < 1, 
Qprcd ~ 1, 



(16) 



n=0 

1.68. ..1.74 for 



with the coefficients {a, 6}„ listed in Table |5| and again the 
uncertainty ctprcd — 1.68... 1.74 emerging from the degener- 
acy of Qprcd = 1 that was mentioned earlieiQ. 

In Fig.|9l the correlation between Qobs/prcd and the ob- 
served low-mass MF slope aobs is assessed for the subsample 
of clusters from the pre sent study that is also considered in 
iDe Marchi eTal] (|2007t) . For comparison, the relation for 
the predicted low-mass MF slope aprcd as a function of 
Qobs/pred IS included as well. Most of the observed data 



Consequently, it represents the same uncertainty as 
*/*dir' = 0...0.2, with Qfprod = 1.68 corresponding to t/td"^' = 
0.2 and aprcd = 1.74 to tjt^l^^^ = 0. Since for most GCs under 
consideration the elapsed fractions of the total disruption time 
are closer to t/^'is"' = 0.2 than to f/td"*""' = (see TableU), we 

adopt Oprcd = 1.68 if Qprcd = 1. 
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Table 5. Coefficients for the fourth-order powerlaw ap- 
proximation of ttpred as a function of Qprcd (see Eq. [T| 
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Fig. 9. Correlation between the observed slope of the low- 
mass stellar mass function aobs and the relative mass-to- 
light ratio with respect to the canonical value Q. The dotted 
curves indicate the theoretically predicted relation between 
a and Q (not a fit), with the dot at the right-hand tip repre- 
senting the canonical values of ag = 1-74 and Q = 1. Values 
of a representing the global MF are marked with triangles, 
while those for clusters with worse than Icr agreement be- 
tween Qobs and Qpred are denoted by squares. Left: For the 
observed mass-to-light ratio fraction Qobs. Right: For the 
predicted mass-to- light ratio fraction Qprcd- 



match the predicted relation between a and Q within their 
error bars, albeit with substantial scatter. This is due to the 
large uncertainties of the observations and possibly also re- 
lated to biases introduced by comparing central and global 
mass-to-light ratios (see Sect. 14. 3|) . The poor quality of the 
observations is illustrated by this spread and by the large 
error bars. For the predicted mass-to-light ratios the trend 
is more well-defined, but for low values of ctobs it does not 
extend down to the mass-to-light ratios that are predicted 
by theory. This could imply that either ttohs or Qpred are 
biased. If the latter is true, it suggests that some GCs per- 
haps dissolve more rapidly than presently included in the 
models. 

As shown in Sect. 14.31 comparison of Ofobs with the ob- 
served and predicted fractions of the canonical M/ Ly ratios 
Qobs/prcd for the GCs with agreement parameter > 2 (see 
Table 3]) provides an independent check of our predicted 
M/Lv ratios. While for NGC 104 and 5272 this does not 
allow for any definitive conclusions, for NGC 6341 and 6809 
the observed mass functions are clearly more consistent 
with our predicted M/Ly ratios than with the observed 
values. 

With Eq. [ini we can also predict the slope of the 
low-ma ss MF fo r clust ers that where not considered by 
De Mar chi et al.l (|2007D . The predicted slopes are hsted in 
Table jSj As stated above, for most clusters with observed 
values of a, the agreement between observed and predicted 
a is reasonable. Only for NGC 6218 and 6712 there is a 
strong discrepancy. For NGC 6218, we expect the devia- 



J. M. D. Kruijssen and S. Mieske: On the reduced mass-to-light ratios of Galactic globular clusters 



13 



MF slopes 



NGC 


Qobs 


<^prcd 


104 


1.2 ±0.3 




288 


0.0 ±0.3 


0.961^3? 


1851 




t-08_(, QQ 


1904 




l.UO_Q 


4147 




1 40+0" 


4590 




1.661H5 


5139 


1.2 ±0.3 


1.55tg-l? 


5272 


1.3 ±0.3 


l-D8_o.oo 


5466 




-, ry+O.OQ 
'-■'Jl -0.12 


5904 






6093 




, r9-|-0.14 


6121 


1.0 ± 0.2 


0.62^°;^3^ 


6171 




i-34l°;l° 


6205 






6218 


-0.1 ± 0.2 


1.44«:?? 


6254 


1.1 ±0.3 


1 50+0.05 


6341 


1.5 ±0.3 


-I qr-l-0.12 


6362 




0.85«i? 


6656 


1.4 ±0.2 




6712 


0.9 ± 0.2 




6779 




'U_o 37 


6809 


1.3 ±0.3 


1 44+004 

i-44_o 06 


6934 




1 fi«+ooo 


7089 




1 fiS+000 
l-D»_o.oo 



Table 6. Observed low-mass stellar niass function 
(MF) slopes ctob s in th e range m = 0.3 — 0.8 M0 from 
iDe Marchi et al.1 (|2007f ) (second column), and predicted 
values Qfpred for the same mass range based on the fit of 
Eq. [Tni and Fig. [5] (third column) . The observed low-mass 
MF slope of clusters representing global MFs from multi- 
mass Michie-King models are denoted in boldface, while 
the other values designate local values close to the half- 
mass radius (see text). The standard errors on aobs is are 
(Ta = 0.2 for the global MFs and ctq = 0.3 for the other 
values (De Marchi, private communication). 



tion to arise from the observed value of aobs, since the pre- 
dicted and observed M/Ly are in excellent agreement (see 
Tabled]). On the other hand, for NGC 6712 the incompati- 
bility may be due to a slight overestimation of (M / Ly)pj-cd 
and thus of Qprcd and Ofpicd- 

In this context it must be noted t h at th e compila- 
tion of a values from iDe Marchi et al.l (j2007t ) is drawn 
from a sample of literature estimates, most based on HST 
data, observed in somewhat different radial regions of each 
cluster. Four of the eleven GCs that coincide with our 
sample of 24 GCs do have a direct estimate for their 
global mass func tion (see Table [Ht - For th e remaining seven 
other GCs from iDe Marchi et al.l (|2007T ). that estimate is 
taken from measurements restricted to the region around 
the half-mass radius rh, of which it is known that the 
shape of t he MF is co mparable to the glo bal (i.e., cluster- 
wide) MF (iRicher et al. 1991: Baumgard t fc Makinall2003t 
iDe Marchi fc Pulond 12007 ) . However, the uncertainty of 
these slopes is larger, and they do not provide a self- 
consistent way to derive the global MF. It is clear that a 
direct determination of the global stellar mass function for 
most of the 24 GCs investigated in this study would allow 
to verify the predictions of the present study with much 
higher confidence. 



In Fig. [TO] we investigate how feasible it is to observa- 
tionally verify the predicted drop of a for our full sample of 
24 GCs. We plot the apparent y-band magnitude of stars 
with 0.3 solar masses Vb.3 for each cluster versus the pre- 
dicted slope apred of the stellar mass function for 0.3 < 
m/MQ < 0.8. Vb.3 is obtained from the distance modulus 
of each GC and from the assumption that Mvo 3 = 9.8 mag 
(jBaraffe et al..,1997i) . We also plot the angular size of each 
cluster versus apred- As a consequence of their generally 
larger galactocentric distance, those GCs with the faintest 
Vb.3 > 25 mag would not be expected to exhibit a strong 
low-mass star depletion. The angular half-mass diameters 
of the GCs with Vb.3 < 25 range between 2 and 8 arcmin- 
utes. To obtain a representative estimate of the global mass 
function, it is clear that wide-field ground-based imaging is 
required for most GCs. For this wide- field imaging, a com- 
pleteness magnitude of V ~ 26 mag is desirable, which will 
allow moderately precise photometry already for y ~ 25 
mag. For 8m class telescopes and with optical seeing in the 
range 0.8 to 1.0", this requires 1-2 hours integration time 
per filter, or 2-4 hours for a two-band exposure. With wide- 
field imagers such as VIMOSOVLT, IMACS@Magellan, or 
SuprimeCam@SUBARU, single-shot images will be suffi- 
cient to cover at least 2-3 half-light radii for most clusters. 
From Fig. [TO] we conclude that the best candidate that also 
complements the compilation by De Marchi et al. (2007) is 
NGC 6779, followed by NGC 6362 and possibly NGC 6171. 

6. Discussion and conclusions 

In this section, we provide a summary and a discussion of 
our results. We consider the effects of the assumptions that 
were made and reflect on the implications of the results. 

6.1. Summary 

In this study, we have investigated the dynamical mass-to- 
light ratios of 24 Galactic globular clusters. We have tested 
the hypothesis of the preferential loss of low-mass stars as 
the main explanation for the fact that the average observed 
mass-to-light ratios of the Galactic GCs in our entire sam- 
ple are only 741^% of the expectations from stellar popula- 
tion models. Accounting for the orbital parameters we de- 
rived dissolution timescales due to two-body relaxation and 
disc shocking for our globular cluster sample and calculated 
the evolution of their masses and p hotometry using the 
SPACE analytical cluster models from iKruiissen fc LamersI 
(2008, throughout this paper KL08). These models account 
for the preferential loss of lo w-mass stars which is fi t ted to 
the A^-body simulations by iBaumgardt fc Makinol (|2003f ) 
and therefore provide non-canonical M/Ly predictions. We 
find the derived dissolution timescales to be in good agree- 
ment with the range required for low-mass star depletion 
to exp lain the observed M/L ratio decrease from lKruiissenl 
(l200l . 

The present-day {t = 12 Gyr) M/Ly ratios 
have been compared to t he ob served values from 
[McLaughlin fc van der Marell (|2005f ). yielding Icr agree- 
ment for 12 out of 24 GCs. We considered possible causes 
for the remaining > Icr discrepancies that occur for the 
other GCs. It is found that 11 of these clusters have pre- 
dicted MjLy very close to the canonically expected M/Ly 
ratios due to their long dissolution timescales and the cor- 
respondingly modest low-mass star depletion, while their 
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observed M/Ly are lower. This is probably due to the 
method by which the observed M/Lv are derived, which is 
biased towards the central M/Ly while our models predict 
global MjLy. For mass-segregated GCs with long disso- 
lution timescales, both values can be substantially differ- 
ent. The discrepant GCs have higher than average King 
parameters Wq, which should inde ed reach mass segrega- 
tion on shorter timescales (see e.g.. iBaumgardt fc Makinj 
l2003f ). This explanation for the discrepancy between some 
of the observed and predicted MjL ratios is confirm ed by a 
study of low-mass star depletion in GCs by JDc Marc hi et al.l 
(|2007l ). whose observed mass functions are in good agree- 
ment with our predictions. The average observed M/Ly 
ratio of 741^% of the canonical expectations would there- 
fore be underestimated. Excluding GCs which likely have 
dissimilar global and central M / L ratios by making cuts 
in dissolution timescale and King parameter, we find that 
the observed and predicted M/Ly ratios are consistent at 
78^11% and 78 ± 2% of the canonical values, respectively. 
For the entire sample, the average predicted fraction of the 
canonical M/Ly ratio is 85 ± 1%. 

To assess the imprint of low-mass star depletion on 
the slope of the low-mass stellar mass function, we 
compared the observe d mass function slopes Oiohs from 
iDe Marchi et all (|2007f ) for 11 GCs contained in our study 
to the values predicted by our models as well as to the 
observed and predicted mass-to-light ratio fractions of the 
canonical values Qobs and Qprod- Most of the measured 
slopes agree with the predictions, but exhibit considerable 
scatter. Since most of them are values derived at around 
the half-mass radius and are extrapolated to global values, 
we also discuss the feasibility of observations for directly 



measuring global mass functions of most of the GCs in- 
vestigated. We show that deep (ground-based) wide-field 
imaging would be necessary, with point source detection 
limits V '^26 mag. The most suitable candidate for such a 
campaign would be NGC 6779. 

6.2. Propagation of assumptions 

In the course of the study presented in this paper, several 
assumptions were made that affect the results to different 
extents. Their implications are as follows. 

(1) We have adopted the SPACE cluster models (KL08), of 
which the stellar evolution and photometry are based on 
the Padova 1999 isochrones (see Sect. [2]). Consequently, 
the predicted cluster photometry and corresponding 
mass-to-light ratios are affected by that choice. To in- 
dicate t he level of the devia t ion w ith the cluster mod- 
els from iBruzual fc CharlotI (|2003D . in Fig. [T] we com- 
par ed the canonically e x pecte d M/Ly from SPACE to 
the IBruzual fc Chariot (|2003f) values for our cluster 
sample. The difference between both is inadequate to 
explain any systematic tendency of low mass-to-light 
ratio with respect to the SPACE models. Therefore, we 
conclude that the adopted cluster models do not effect 
substantial implications for the predicted M/ L ratio^. 

(2) Due to the treatment of the preferential loss of low-mass 
stars in the SPACE cluster models, there are indications 



Of course, if the preferential loss of low-mass stars is ac- 
counted for, i.e., non-canonical models are considered, the SPACE 
cluster models predict very differ ent photometric evo l ution than 
canonical cluster models such as IBruzual fc Charloti (|2003f l. 



J. M. D. Kruijssen and S. Mieske: On the reduced mass-to-light ratios of Galactic globular clusters 



15 



that the predicted mass-to-light ratios could be underes- 
timated during the final ^ 15% of the total cluster life- 
time (KL08). Table m shows that none of the GCs in our 
sample reside in this regime. On the other hand, we did 
not include primordial mass segregation, which could 
decrease the predicted M/L ratios by ~ 10% (KL08). 

(3) By adopting the average cluster orbits from 
iDinescu et alT ()1999D . we assume constant orbital 
parameters over the total cluster lifetimes. Considering 
the ballistic nature of the orbits, such an assumption 
is legitimate as long as the external conditions do 
not strongly differ. The Galactic potential was only 
substantially different from its present state during 
the formation of the Milky Way. A more extended 
distribution of mass during these early epoch would 
obviously increase the dissolution timescale due to disc 
shocking, and would affect the dissolution timescale 
due to two-body relaxation in a similar way because 
of the reduced tidal field. Consequently, this would 
imply that the mass loss during the first ~ 1 Gyr of 
our models is overestimated, causing our initial masses 
to be overestimated as well. However, the extended 
nature of the Milky Way would cause dissolution 
due to giant molecular cloud en counters to becorn e 
an important mechanism (e.g., iGieles et al.l l2006f) . 
thereby counteracting the previous effect. Although we 
cannot rule out any consequences, a residual influence 
would only be relevant for a small fraction (the first 
~ 10%) of the total cluster lifetime, where mass loss by 
dissolution is much less effective than later on during 
cluster evolution. Therefore, this likely only affects our 
analysis within the error margins. 

(4) We have compared our predictions to the observed 
mass- to-light ratios from iMcLaughlin fc van der Mare] 
(2005), which are biased towards central M/L values. 
As treated more extensively in Sect. 14. 3[ this yields un- 
derestimated observed M/L ratios for mass-segregated 
clusters with long dissolution timescales. Therefore, 
based on the earlier discussion and the parameter range 
in which discrepancie s arise, we consider the dynamica l 
(M/iy)obs from .McLaughlin k, van der Marell (|2005[) 
to be subject to improvement for GCs with both dis- 
solution timescales tot > 5 Myr and King parameters 
Wq>7. 

6.3. Consequences and conclusions 

The consequences of our findings are not only relevant to 
studies of the M/L ratios of compact stellar systems, but 
also to other properties of these structures. Here we list 
them together with the conclusions of this work. 

(1) When constraining our sample to the subset for which 
the observed M/Ly likely reflect the global values, we 
find that the preferential loss of low-mass stars can ac- 
count for the ^20% discrepancy between observed dy- 
namical mass-to-light ratios of Galactic GCs and those 
expected from stellar population mode ls that assurn e 
a canonical present day mass function (|Kroupall200lh . 
This alleviates the factor of two offset in M/L between 
GCs and UCDs by about 25%. Still, some additional 
dark mass with respect to a canonical IMF is required 
to explain the M/L of most UCDs. 



(2) Accounting for the orbital parameters, present-day 
masses and chemical compositions of individual clus- 
ters, we find that there is good agreement between our 
model predictions and observations of the M/Ly ra- 
tios of these clusters. For the GCs with worse than Itr 
agreement there are strong indications that the discrep- 
ancy is due to an underestimation of the observed M / L 
ratio. In mass-segregated clusters with long dissolution 
timescales, the observed M/Ly ratios represent central 
values that do not reflect the global M/Ly ratio. 

(3) The ideal way to confirm the validity of our explana- 
tion for the reduced M / L ratios of GCs for individ- 
ual clusters will be to obtain a homegeneous set of 
deep wide-field imaging for most GCs. This would ex- 
pand and complement the currently available heteroge- 
neous data sets of space based GC imaging, which is 
restricted to small fields in each GC, at different radial 
ranges. By this, the global mass and luminosity func- 
tions could be measured directly for individual GCs and 
be compared quantitatively to the predictions of this 
paper regarding the low-mass star depletion due to dy- 
namical evolution. In addition, velocity dispersion mea- 
surements would allow for the determination of global 
M/L ratios, thus providing an up date to those from 
IMcLaughlin fc van der Marell (|2005f ). 

(4) The topic of globular cluster self-enrichment and mul- 
tiple stellar populations can also be considered within 
the framework of this paper. In a recent study by 
iMarino et ^ (|2008D it is shown that NGC 6121 con- 
tains two stellar populations that are probably due to 
primordial variations in their respective chemical com- 
positions. It is mentioned that the present-day mass of 
NGC 6121 is an order of magnitude smaller than that 
of known multiple-population GCs such as NGC 1851, 
2808 and 5139. Conseguentlv. lMarino et all 12008') pose 
the question how the enriched material could have re- 
mained in such a shallow potential and argue that mul- 
tiple populations are unlikely to be strictly internal to 
GCs, unless they are the remnant of much larger struc- 
tures. In the case of NGC 6121, our calculations seem 
to explain the issue, as it is the initially fifth most mas- 
sive GC of our sample (Mdj 10^ Mq). As a result, 
mass could have been retained much more easily, imply- 
ing that the multiple populations of NGC 6121 are no 
reason to invoke external processes for enrichment and 
to abandon the self-enrichment scenario. 

We conclude that the variation of M / L ratio due to clus- 
ter dissolution and low-mass star depletion is statistically 
significant and serves as a plausible explanation for the 
difference between observed and canonical M/L ratios. 
Moreover, it has several implications that should be ac- 
counted for in GC studies, since its effects can be accu- 
rately quantified. We also suggest that the M/L decrease 
is considered in independent observational verifications to 
further constrain the evolution of the stellar mass function 
in dissolving globular clusters. 
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Appendix A: Error analysis 

In this Appendix, the error propagation through our com- 
putations is discussed. The errors in Tables [2] and [3] are 
standard errors, most of them determined by computing 
the formal error propagation. For a function f{xi, X2, Xi) 
this implies 

with (Ti the error in the parameter i. Asymmetric errors on 
each parameter are both separately propagated by employ- 
ing the same recipe, while inverse relations are accounted 
for by swapping the positive and negative errors. However, 
Eg. I A. II assumes an approximately constant derivative over 
the standard error interval. For very large errors on non- 
linear relations this assumption does not hold. The first of 
two parameters where we have to correct for this effect is 
to fih ■ It is inversely r elated to the destruction rate Vgh from 
iDinescu et al.l (|l999t) . which is a parameter with very large 
relative errors, even to the extend that after computing the 
error propagation one can have io,sh — <^to sh ^ ^- Because 
dissolution timescales smaller than zero are not physical, 
instead the negative error on io,sh is determined by com- 
puting 

'^to.sh ^ ^0,sh(l^sh,M"cl,a:;corr) (A. 2) 

where a~ indicates the negative error in a parameter i. In 
the context of Eg. lA.li this approach is equivalent to assum- 
ing the derivative equals the mean slope of f{x) over the 
interval [x — Ca;, x]. On two occasions (NGC 6093 and 6712), 
a strongly asymmetric error in ^o.sh propagates into to. tot 
such that io,tot+crttt„t > min(to,cvap+cri+ _^^^p,^o,sh+o•fl; ,J• 
However, since a very large positive error in io,sh or io,cvap 
would make the term vanish in the inverse addition of Eq.[Bl 
it should not propagate into a similarly large error in to, tot- 
This brings up the second parameter we have to correct for 
the propagation of large errors through non-linear relations. 
We define the error in to, tot for NGC 6093 and 6712 such 

that to,tot + Crto,,^t = *0,evap + CTtt,,„,p • 

The error margins on our predictions in Table |4] are 
determined by numerically evaluating Eq. lA.ll for the de- 
sired quantities. Our predictions depend on the observed 
mass, metallicity and dissolution timescale. The derivatives 
of M I Ly with respect to the former two are trivial since 
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M/Lv is determined by interpolating over these parame- 
ters. For the dissolution timescale, we compute additional 
models at tg.tot — cr^" to obtain the numerical deriva- 
tive of MjLy with respect to io.tot- In fact, this is the 
differential rather than the derivative, because for long dis- 
solution timescales MjLy can be locally constant, while 
it varies over a larger range. The only case were a non- 
linearity forces us to derive alternative errors is for the 
positive standard error on M/Ly- Although the uncer- 
tainty in metallicity could increase the predicted mass-to- 
light above its canonical value, the uncertainty in mass and 
dissolution timescale cannot due to the flattening of the 
cluster isochrones in the {M, M/Ly}-plane (see Fig. U]). 
Therefore, the combined positive standard error of the 
mass-to-light ratio due to the uncertainty in mass and dis- 
solution timescale crt;^'*" is defined as 



+ ,M,-to 



M/Lv 



Ki/L.^iM/LyU^ - {M/Ly\,,^ , (A.3) 



with o"!^'*^'^" the standard error according to Eq. lA.ll 

(M/iy)can thc canouically expected mass-to-light ratio 
and {M/Ly)p^cd the predicted value. This definition en- 
sures that the positive standard error is never larger than 
the difference between the canonical and predicted mass- 
to-light ratios. 

Except for the alternative error in Eq. IA.3I that is spe- 
cific to M/Ly, the standard errors on the predicted initial 
masses are determined analogously to the above. For the 
remaining lifetimes, numerical derivatives with respect to 
mass, metallicity and dissolution timescale are simply ob- 
tained by reintegrating Eq. [T] for slightly different initial 
conditions. 

Finally, for the predicted slopes in TablelHl the errors are 
computed using Eq. lA.ll and restricted such that a + ctq < 
ao (analogous to Eq. IA.3|) . 



